Periodicity of the time-dependent Kohn-Sham equation and the Floquet theorem 
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The Floquet theorem allows to reformulate periodic time-dependent problems such as the inter- 
action of a many-body system with a laser field in terms of time-independent, field-dressed states, 
also known as Floquet states. If this was possible for density functional theory as well, one could 
reduce in such cases time-dependent density functional theory to a time-independent Floquet den- 
sity functional theory. We analyze under which conditions the Floquet theorem is applicable in a 
density-functional framework. By employing numerical ab initio solutions of the interacting time- 
dependent Schrodinger equation with time-periodic external potentials we show that the exact effec- 
tive potential in the corresponding Kohn-Sham equation is not unconditionally periodic. Whenever 
several Floquet states in the interacting system are involved in a physical process the corresponding 
Hartree-exchange-correlation potential is not periodic with the external frequency only. Using an 
analytically solvable example we demonstrate that, in general, the periodicity of the time-dependent 
Kohn-Sham Hamiltonian cannot be restored by choosing a different initial state. Only if the external 
periodic potential is sufficiently weak such that the initial state of the interacting system evolves 
adiabatically to a single, field-dressed state, the resulting Kohn-Sham system admits the application 
of the Floquet theorem. 

PACS numbers: 31.15.ee,31.15.ec,42.50.Hz,32.80.Rm 



I. INTRODUCTION 



The properties of multi-electron systems can in prin- 
ciple be predicted by solving the interacting many-body 
Schrodinger equation. However, numerical solutions are 
only feasible for small systems consisting of a few in- 
teracting electrons due to the exponential scaling of the 
computational demand with the number of particles. 
One possible way to overcome this so-called "exponential 
wall" [1^ is density functional theory (DFT) which 
has been successfully applied to many-body systems in 
a wide range of areas in physics and chemistry. DFT 
is based on the existence of an energy functional whose 
minimization yields the ground state density. This min- 
imization is usually performed via the so-called Kohn- 
Sham (KS) construction [4 where the interacting multi- 
particle system is mapped to a unique system of non- 
interacting particles having the same ground-state den- 
sity. The non-interacting problem decouples into (non- 
linear) one-particle equations with an effective Hamilto- 
nian depending on the density. The advantage of DFT 
thus originates from the fact that the solution of N one- 
particle equations is less involved than solving one ex- 
ponentially scaling A'^-particle problem. The crucial in- 
gredient in the KS construction is the Hartree-exchange- 
correlation (Hxc) potential which, if it was exact, would 
comprise all many-body effects. Its exact form is, how- 
ever, unknown in general, and approximations have to be 
employed in practice. 



Time-dependent density functional theory (TDDFT) 
extends DFT to time-dependent problems [Ml] • The ex- 
istence of a time-dependent KS system is, however, no 
longer based on a minimization principle but on the lo- 
cal force equation of quantum mechanics E] • The ex- 
act time-dependent Hxc potential has additional, subtle 
features: it depends on the initial state (both interact- 
ing and non-interacting) and on the density at previous 
times (that is, it has "memory") j^]. Therefore the con- 
struction of better approximations to the Hxc potential 
within TDDFT is much more involved than in DFT. This 
is even more of a problem as it turns out that the Hxc 
approximations known from DFT often fail when applied 
to TDDFT beyond linear response |10|. 

One might think that if TDDFT is employed to the 
study of multi-particle systems subject to time-periodic 
external potentials, e.g., an atom interacting with a 
monochromatic laser field, one could involve the Floquet 
theorem. Indeed, in such situations the interacting many- 
body time-dependent Schrodinger equation is a partial 
differential equation with time-periodic coefficients and 
thus admits a time-periodic basis. As a consequence, the 
problem can be converted into an infinite set of time- 
independent equations by virtue of the Floquet theorem 

Already at the beginning of the application of density- 
functional theory to time-dependent systems attempts 
were made to incorporate Floquet theory in a density- 
functional framework [15j . A minimization principle was 
proposed, which was perturbative in nature and hence 
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valid only for weak and ofF-resonant fields. However, 
even if these conditions are met there are problems with 
defining a proper adiabatic limit, which is fundamental 
to the proposed minimization procedure |16j . The prob- 
lems arise due to the fact that Floquet theory maps the 
quasi-spectrum of the time-dependent problem into an 
interval of length w, i.e., the frequency of the periodicity 
employed. In any interval I = {x — uj/2, x -|- w/2} for x e 
M arbitrary we find infinitely many quasi-eigenenergies 
(they are dense in /) and thus infinitely many eigenfunc- 
tions around every point in the quasi-spectrum. A con- 
sequence of this is that there is no unique final state to 
which the system tends as the external perturbation is 
turned off adiabatically. In order to restore the adiabatic 
limit a truncation to a finite basis is usually employed, 
which is anyway unavoidable in practical calculations. 

In Refs. p!7H20] Floquet-DFT approaches were pur- 
sued for non-perturbative fields and later criticized in 
Refs. [1T1[31] where the authors also suggested to embark 
upon the problem from a TDDFT point of view, thereby 
avoiding the minimization problem. The basic question 
then remains whether a Floquet basis can be found for 
the associated KS system, i.e., whether the KS Hamilto- 
nian itself is periodic. Known explicit expressions for the 
exchange-correlation potential in the time-dependent KS 
Hamiltonian such as the adiabatic local density approxi- 
mation or generalized gradient approximations [3J |B] have 
the feature that a periodic density will lead to a periodic 
KS Hamiltonian (with the same period) since the adia- 
batic Hxc potentials depend on the instantaneous density 
only. However, the density does not have to be periodic 
and, in fact, it generally is not, as we will demonstrate 
in this work. On the other hand, even if an approximate 
functional leads to an aperiodic KS potential because of, 
e.g., an aperiodic density, this does not yet demonstrate 
the incompatibility of TDDFT and Floquet theory, be- 
cause the unknown exact KS potential nevertheless could 
be periodic. In this work we will show by means of nu- 
merical and analytical counter examples that this, unfor- 
tunately, is not the case and thus TDDFT is, in general, 
not compatible with Floquet theory. 

The paper is structured as follows. In Sec. [IT] we review 
the basics of Floquet theory from a TDDFT perspective. 



In Sec. HI we compute the exact KS potential for a two- 
electron model system, present the Fourier-transformed 
exact KS potential, and investigate whether the Floquet 
theorem is applicable to the KS Hamiltonian. In Sec. IV 
an analytical example is given to analyze the initial-state 
dependence of the KS potential and its relation to the pe- 
riodicity of the KS Hamiltonian. We conclude in Sec. |V] 
For simplicity, we restrict ourselves to one-dimensional 
systems in this work. Such systems are frequently used 
in the theory of laser-matter interaction because they 
can be solved numerically exactly, and they are known 
to capture many of the essential features of their three- 
dimensional analogs. All equations in this work can 
be straightforwardly extended to the three-dimensional 
case. 



Atomic units h = me = |e| = Aire^ = 1 are used 
throughout unless stated otherwise. 



II. BASIC THEORY 

Consider a system of N interacting electrons governed 
by the Hamiltonian 



H{t) =f + V,, + V{t) 



(1) 



with, in position-space representation, the kinetic energy 
operator 



N 



^ 2dxl 



(2) 



the interaction potential 



-y 

2 ^ 

and the external potential 



^nn — 2 2^ Vcc(\Xi — Xjy, 



(3) 



(4) 



We assume the interaction to be Coulombic. In one- 
dimensional models the Coulomb-interaction is usually 
smoothed by a softening parameter e > 0, 



'^cc I ) 



1 



^/{Xi - Xj)^ + < 



(5) 



We further specialize on external potentials consisting of 
the interaction with a (static) nucleus of charge Z and a 
laser field E(t) in dipole approximation, i.e.. 



v{xi,t) 



(6) 



The eigenstates and eigenenergies of the laser field-free 
system at time t = are obtained via the solution of the 
time-independent Schrodinger equation 

H{t)'i'{xiai ■ ■ ■ xn^n) = E'^{;xiai ■ ■ ■ xn^n)- (7) 

Here '^{xi<ti - ■ ■ xjqUM) is an antisymmetric A^-particle 
eigenfunction of the space and spin variables , Ui, and £ 
is its eigenenergy. In order to obtain '^(xiUi • ■ ■ xn<Jn, t) 
for t > one may solve the time-dependent Schrodinger 
equation (TDSE) 

idt'^ixiCTi ■ ■ ■ XNcrN,t) ^ i/(<)*(a;i(Ti • • • XN<yN,t) (8) 

for a fixed initial state 4'o(a;icri • • • XNa^). However, due 
to the "exponential wall" [T] it is computationally very 
challenging to solve this equation. In fact, in the case of 
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intense laser fields where the numerical grids need to be 
large it is feasible only for A'' < 3. 

Now we turn our attention to the non-interacting KS 
system that, by construction, yields the same single- 
particle density n(x, t) as the interacting system. For 
simplicity, we assume that we are dealing with spin- 
neutral systems. The KS Hamiltonian then reads 



with 



n{t)^HKsi[n];t)~idt 



(16) 



HKs{[n];t) = 



1 



i.e., ^ct assumes the role of an eigenvalue and 4>aix,t) is 
the corresponding eigenstate. If so, also 

(17) 



-v{x,t) -f whxc (["■]; a:, t), (9) are solutions of the eigenvalue equation (15). Owing to 



where v{x,t) is the external potential ^ and 
VY{-Kc{[T^]]X,t) is the Hxc potential which is a functional 
of the single-particle density n{x,t) (for notational sim- 
plicity we do not indicate the dependence on the initial 
states) . The two potential terms combined are called the 
KS potential, i.e., 



VKsi[nhx,t) = v{x,t) +vii^ci[n];x,t). 



(10) 



In what follows we assume that the external laser field 
is monochromatic with a period wi . The time-dependent 
KS equation reads 



idt^k{x,t) = HKs{[n];t)'S?k{x,t), 



(11) 



where ^k{x,t) is the k-th KS orbital for the KS parti- 
cle with initial state ^k{x,0)- The time-dependent one- 
particle density n{x,t) then is 



N 



(12) 



k=l 



Now we make the basic assumption of any Floquet ap- 
proach in a density-functional framework: if the Hamil- 
tonian H{t) describing the N interacting electrons is pe- 
riodic with the frequency wi, i.e., E{t + T) = E(t) with 
T = 27r/a;i, then we assume the same periodicity for the 
KS Hamiltonian as well. We neglect for the moment po- 
tential problems with respect to the non-linear nature of 
the KS equations, which will be discussed in detail in the 
subsequent Sections of this work. 

If the KS Hamiltonian is periodic with T then, by 
virtue of the Floquet theorem, we can write the KS or- 
bitals in a time-periodic (Floquet) basis {(/)q,(x, t)}^^^ as 



<^k{x,t) ^^CkaG '^•'*<l>a{x,t), 



(13) 



where the are the so-called quasi-energies and Cka — 
{4>a{t = 0)\^k{t = 0)). Further, the (j)a{x,t) are periodic 



in T, i.e.. 



(t>a{x,t) = (j>aix,t + T). 



(14) 



The Floquet orbitals (t>a{x,t) fulfill the eigenvalue equa- 
tion 



'H{t)cj)a{x,t) ^^a(l)a{x,t) 



(15) 



the time periodicity of (t)a{x, t) we can write 

0„(a:,<) =^^„,,(a:)e-'''^i*, / G Z. (18) 



With Eqs. ( 13 ) and ( 18 1 the KS orbital can thus be writ- 
ten as, 



y^a,lix), 



(19) 



where the eigenstates {fa.i{x)}aeK,iei. form the time- 
independent Floquet basis. 

We divide the Hamiltonian HKs{[n'];t) into a time- 
independent part 



Ho 
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Z 



(20) 



2dx^ Vx^' 

the coupling to the monochromatic external field 

xE{t) = i;+(a;)e'"i* + v-{x)e-''^'\ (21) 

and V}i-xc{[n\; x,t). Since we tentatively assume time- 
periodicity of the whole KS Hamiltonian we can write 

«Hxc(W;a;t) = ^e-""i*Kxc(N;a;)]i, (22) 



I £ Z. Plugging the expansions (frsl), (21), and (22 1 in 



Eq. ( 15 ) we obtain the TDDFT-Floquet equations ,17 



{^a + 1^1 - Ho)ipa.l{x) (23) 
= V+{x)ipa.l-l{x) + {x)(paJ+l{x) 

+ ^{vH^c{[n];x))i_rn{x)(Pa,l{x). 

m 

The index I of a Floquet state (pa,i{x) is known as the 
"block index," which may be interpreted as the number 
of photons involved in the process under study (provided 
one arranges that the / = 0-block adiabatically connects 



to the field-free situation). The Floquet equation (23) 



couples any Floquet block / to its neighboring blocks / ± 1 
via absorption or emission of a photon. Contributions of 
non-neighboring blocks may only be included through the 
Fourier-components of the Hxc potential. This is differ- 
ent from the Floquet equations for the interacting TDSE 
which couple only neighboring blocks because E{t) is the 
only time-dependent element in the TDSE Hamiltonian. 
However, in the TDSE case the Floquet basis functions 
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depend on all spatial variables, not just on a single one 
as in the KS case. 



In principle, Eq. (23) is an infinite-dimensional set of 



coupled partial differential equations, in practice, it is 
truncated so that /min < I < ^max where |^min| and |Zmax| 
should be large enough to capture all the relevant pro- 
cesses in which photons are emitted or absorbed. 

If Eq. ( 23 ) was valid, the periodic time-dependent 



many-body problem would be significantly simplified be- 
cause the time-dependence had been eliminated via Flo- 
quet theory and the "exponential wall" via DFT. 



III. PERIODIC OR APERIODIC KS 
HAMILTONIAN? 

In order to prove that Floquet theory is generally not 
applicable to TDDFT it certainly is sufficient to find one 
counterexample. However, a Floquet approach might 
still be useful as an approximative approach, especially 
given the fact that TDDFT in practice is itself approxi- 
mative anyway. Hence, we analyze under which circum- 
stances the KS Hamiltonian is periodic or not. In or- 
der to do so we employ a widely used numerically ex- 
actly solvable one-dimensional model Helium atom [23l - 
[25] . In this model both electrons move along the laser- 
polarization direction only, and the Coulomb interac- 
tion is replaced by a soft-core potential as introduced 
in Sec. |llj The TDSE Hamiltonian of the model system 
thus corresponds to the Hamiltonian Q with N = 2. 
The smoothing parameter was e — 1, as, e.g., in [25] . 

The initial TDSE state is chosen to be the spin-singlet 
ground state of the interacting system 

^'o(a;iai,a;2CT2) = ■^o{xi,X2)-^ (| ti>l i2) - I W) ■ 

(24) 

Since the Hamiltonian is spin-independent, the system 
remains also during the dipole interaction with a laser 
field in a spin-singlet configuration, and we can concen- 
trate on the symmetric spatial part "^q^xi, X2) of the 
wave function only. The TDSE ([8| is solved numeri- 
cally using the Crank-Nicolson propagator to obtain the 
time-dependent spatial wavefunction '^{xi,X2-,t). 

In Ref. [Tl] we introduced a method to extract the pop- 
ulated Floquet states of the interacting system directly 
from '^{xi,X2,t). By controlling the laser parameters 
we can either have an adiabatic evolution of the field- 
free state 'i'o{xi,X2) to a field-dressed (Floquet) state 
or a non-adiabatic one, where several Floquet states are 
populated. The laser intensity, frequency and the ramp- 
ing time decide on the adiabaticity of the time-evolution 
of the interacting system. For adiabatic evolution we 
have in the TDSE-Floquet calculation only one relevant 



Floquet-state index a in the TDSE analog of (19), 



«'(xi,X2,<) = I^c„e-'«"*^e-""iVa.i(a^i,a;2). (25) 



Hence, in this case 



^'(a;i,X2,i)~e"'«°*^e-''"^* 



^a,lixi,X2), 



(26) 



and the density n{x,t) = 2/ dx' \^ {x , x' , t)]"^ will only 
have frequency components proportional to multiples of 
the laser frequency uji . The KS Hamiltonian depends on 
the density. If the KS potential is periodic with respect to 
integer multiples of the laser frequency there would be no 
problem because HKs{[n{t + T)];t + T) = HKs{[n{t)];t), 
and thus the Floquet theorem still holds. Instead, frac- 
tional harmonics or, even worse, incommensurate fre- 
quencies in HKs{[n'{t)];t) would render the Floquet the- 
orem inapplicable. If more than one Floquet state is 
populated, say a — ai and 02, the Fourier-transformed 
density n{x,uj) will also have frequency components pro- 
portional to the inverse of the quasi energy difference 
ICaa ~ CqiI- It would be mind-bogghng if the unknown 
exact Wxc(["-] ; x) was able to remove such frequencies from 
H-i^g([n{t)];t). However, in order to prove that in gen- 
eral the exact Wxc ( [n] ; x) contains frequency components 
different from uji we construct the exact Uxc ( [n] x) ex- 
plicitly in the following for both the adiabatic as well as 
the non-adiabatic evolution of the field-free state to the 
field-dressed states. 

Once we have obtained 5'(a;i, a;2, i) by solving the 
TDSE ([s]) we can construct the exact KS orbital and 
the potential following Refs. [Ml US]- In the two-electron 
spin-singlet case the KS wave function consists of only 
one spatial orbital $(x,t), i.e., 

^{xiai,X2(T2,t) 

= Hx,,mX2,t)^{\tl)\i2)-\U)\h))- 

The KS orbital can be written as 



(27) 



where n{x, t) is the exact particle density and S{x, t) is 
the exact phase of the KS orbital. The expression for 
the phase in terms of density is given by the continuity 
equation as [26l [27] 



(28) 



-dx [n{x,t)dxS{x,t)] = dtn{x,t). 



Equation (11) can be inverted to write the KS potential 



in terms of the KS orbital as 



idt'S>{x,t) + ^dl<i>{x,t) 



= l^^^^^~d,S{x,t)-Ud.S{x,t)]\ (29) 
2 Vn{x,t) 2 

The imaginary part of the potential is zero due to the 
continuity equation (28). The density n{x,t) and the 
phase S{x,t) are computed from '^{xi,X2,t) and by 
the above construction we obtain the exact KS potential. 
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Such a straightforward construction is possible only if we 
have a single spatial orbital. In the general case of several 
KS orbitals one would need to employ a computationally 
more demanding fixed-point method, as demonstrated in 
Rcfs. [211 m]. Once the exact KS potential is computed, 
it is Fourier transformed in time to investigate its peri- 
odicity. 

Besides the basic problem of the periodicity of the KS 
potential for a given interacting density, there is the in- 
herent non-linearity of the KS scheme. Even though the 
exact KS potential might be periodic for a certain prob- 
lem, it is far from obvious that one can employ a Floquet- 
based KS scheme to predict it. For instance, although 
an adiabatic approximation, e.g., in the two-electron 
spin-singlet case the exact exchange-only approximation 
^(cxact)^j^j.^^ = Jdx' [nix',t)/2]/^{x- x')^ + e, does 
inherit the periodicity of the density, it is not guaran- 
teed that the non-linear KS equations produce a periodic 
n{x, t). This becomes obvious when we consider the iter- 
ative solution of the KS equations, where we start with 
an initial guess for the density that is periodic with wi. 
We then have a periodic KS Hamiltonian from which we 
can (since in every iterative step we have a linear par- 
tial differential equation) infer a Floquet basis. We then 
solve the resulting linear equations and obtain a new den- 
sity. This density will in general not be periodic and we 
no longer find a Floquet basis with period uji only. This 
makes the problem of the non-linearity in connection with 
a Floquet approach evident. 



A. Adiabatic and periodic example 

First we consider an 800-nm (wi — 0.056) laser pulse 
with two cycles ramp-up and 16 cycles of constant am- 
plitude. The electric field amplitude is E = 0.063, corre- 
sponding to a laser intensity of 1.4 x 10^^ W/crn^. It 
turns out that in this case the density dynamics are 
periodic with the laser period. In Fig. [l] we plot the 
exact |t;Ks(2;,w)p over four orders of magnitude vs the 
harmonic order w/wi. Only harmonics of the laser fre- 
quency at all space points of the KS potential are visi- 
ble. The Floquet theorem is applicable in this case, as 
HKs{[n{t + T)];t + T) = H]^s{[n{t)];t) to a high degree 
of accuracy. 



B. Non-adiabatic and aperiodic example 
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FIG. 1. (Color online) Logarithmically scaled plot of 
|i'Ks(2;,t<j)P for u>i = 0.056, E — 0.063, two-cycle ramp-up, 
and 16 cycles constant amplitude. Notice that only harmon- 
ics of the laser frequency are present over a dynamic range of 
four orders of magnitude. Superpositions of Floquet states do 
not play a role, the dynamics are sufficiently adiabatic, the 
Floquet theorem is applicable to H-Ks{[n{t)];t). 
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FIG. 2. (Color online) As Fig.[T]but for oji = 2.6, E = 0.34, 
four-cycle ramp-up, and 172 cycles constant amplitude. The 
Fourier-transformed potential displays anharmonic frequency 
components, i.e., it is aperiodic. 



As a second example we chose a short-wavelength 
17.5-nm (loi — 2.6) laser pulse with four cycles ramp- 
up and 172 cycles of constant amplitude. The electric 
field amplitude E ~ 0.34 corresponds to a laser inten- 
sity of 4 X 10^^ W/cm^. The fast ramping induces a 
non-adiabatic time-evolution and results in a superposi- 
tion of Floquet states in the TDSE result. The exact 
KS potential oscillates with periods related to the in- 
verse of the quasi energy differences. In Fig. [2] this new 



timescale manifests itself as side bands around the multi- 
ples of the laser frequency. The quasi energy differences 
are determined by the field-free spectrum of the system 
under study and by the ac Stark shifts so that it may 
well happen that they are irrational fractions or multi- 
ples of oji. In that case even a T' > T = 2tt/lui for which 
HKs{t + T') = i?Ks(0 does not exist. 
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Resonant interaction 



When the laser is tuned to the exact resonance be- 
tween the initial (ground) state and a dipole-accessible 
excited state, Rabi-oscillations set in, typically on a time 
scale that is much longer than the laser period so that 
for the Rabi frequency $7 one has ^ ui. In this 
case the density is periodic with the Rabi- frequency Q, 
not with the laser frequency uji. At time T1/2 = Tr/fi 
the upper state is populated, at time 27r/f2 the initial 
state is populated again. The Rabi-frequency depends 
on the electric field amplitude E of the laser and the 
transition dipole matrix element /ioi between the two 
bound states involved. Rabi-oscillations are not cap- 
tured in TDDFT when known and practicable adiabatic 
exchange-correlation potentials are used. Of course, the 
density dynamics between the two states are correctly de- 
scribed when the exact KS potential is used, for instance 
for the numerically exactly solvable model-He system em- 
ployed in this work. It is known that after the time 
when the single-particle density n{x,Ti^2) is that of the 
excited interacting system, the exact KS potential is the 
ground state potential to that density |251 130| . In fact, 
there is no stationary state in the KS potential to which 
the population may be transferred. Hence, the exact KS 
system governs the dynamics by an "adiabatic deforma- 
tion" of the ground state density. Despite this extremely 
simple "Rabi-flopping" dynamics, resonant interactions 
are among the worst cases for TDDFT with known and 
practicable exchange-correlation potentials. 



It is well known that a Floquet treatment of the TDSE 
leads to avoided crossings of the two field-dressed state 
energies when plotted as a function of laser frequency 
|13j . At exact resonance the two Floquet states are 
equally populated and separated in energy by hfl. Hence, 
resonant interaction is a prime example where a super- 
position of Floquet states plays a role even if the laser 
pulse was turned on adiabatically. 



The laser frequency in our model simulation was tuned 
to be at resonance between the ground spin-singlet state 
and the first excited spin-singlet state of the model He- 
lium atom, Lo = £1 — £q = 0.533 [3S]. For the chosen 
field amplitude E = 0.016 (corresponding to a laser in- 
tensity of 9 X 10"'^^ W/cm^) the ground state population 
reaches zero at T1/2 ~ 174, i.e., = 0.018. Figure [3] 
shows |wks(2;,w)P for two cycles ramp-up and 148 cy- 
cles of constant amplitude. The Fourier-transformed po- 
tential shows strong sideband peaks at qcdi ± pfl with 
q,p G Z. Hence, while in the previous example of non- 
adiabatic ramping one might argue that the anharmonic 
peaks in the spectra are weak and therefore could be ig- 
nored, a resonant interaction generates sideband peaks 
of strengths comparable to the harmonics. 




1.0 1.5 2.0 
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FIG. 3. (Color online) As Fig. [T] but for the resonant in- 
teraction with to = £1 — £0 ~ 0.533, E = 0.016, two cycles 
ramp-up, and 148 cycles constant amplitude. Peaks at posi- 
tions quji ± pil with q,p £ Z are seen. 



IV. INITIAL STATE CHOICE 

For the above examples of non-adiabatic ramping or 
resonant interaction a minimization procedure with a fi- 
nite Floquet basis would lead to a laser-aperiodic KS 
Hamiltonian that renders the Floquet theorem inappli- 
cable in the first place. From the TDDFT perspective 
we obtain a laser-aperiodic KS Hamiltonian because of 
the time evolution starting from the chosen initial state. 
However, in TDDFT HKs{[n{t)]; t) should actually read 

HKs{[n'{t),'i'o,{^ok}]',t) because of the dependence of 
the time-dependent KS potential on both the interacting 
initial state and the KS initial states $o/c [SlIHI- Thus 
a loophole for a most stubborn assumable proponent of 
TDDFT-Floquet theory remains: a different choice of 
initial KS states ^ok{x) = ^k{x,t = 0) could keep the 
KS Hamiltonian periodic in T. In this Section we give a 
counter example for which all possible initial states lead 
to laser-aperiodic KS potentials if the density is laser- 
aperiodic. To do so analytically we construct a KS sys- 
tem of two non-interacting electrons on a quantum ring 
of diameter L so that we have periodic boundary con- 
ditions [27j. This makes it an ideal system to analyze 
the time-periodicity of the KS Hamiltonian for the vari- 
ous possible initial states. For spin-singlet states of these 
electrons one can describe the system by a single KS or- 
bital (27) as in our model Helium system above (in the 
limit L 00 the quantum ring becomes equivalent to 
the Helium model). Following the procedure outlined in 
Sec. |III| the potential can be written in terms of the den- 
sity and the phase of the KS orbital as 



VKsiim, n],x,t) 



1 dl^Jn{x,t) 

2 ^n{x,t) 



(30) 

dtS{[m,n],x,t) - l-[dxS{[m,n], x,t)]'^ , 
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which is an exphcit functional of the density n = n{x, t) 
and an integer number m € Z. As shown in Ref. |27j . for 
periodic boundary conditions the phase S can be written 
in the integral form 



S{[m,n],x,t) ^ dy Kt{x,y)dtn{y,t) 
Jo 

2nm dz 



(31) 



Jo n(z,tY 

Jo n(z,t) ^ ' ^ 

where the Green's function, Kt{x,y) is defined as 

K,ix,y)^l[Oiy-x)^Oix^y)]f£-^ 
_ V{x,t)i]{y,t) 

fL dz ' 
Jo n(z,t) 

with 9 the Heaviside step function and 

dy , r dy 



ri{x,t) 



n{y,t) 



+ 



n{y,t) 



(32) 



(33) 



Since the KS orbital obeys the periodic boundary condi- 
tions, the phase S has to satisfy 



S{L,t) ^S{0,t) + 2Trm 
d^Smt)^ d^S{0,t). 



(34) 
(35) 



Hence the integer m plays the role of labeling all the 
possible KS orbitals (for different initial-state choices) 
that are consistent with the density n(x,t). 
If we assume that 



n{x, t + T) — n{x, t), 



we have 



dz 



X n{z,t) 



dz 



n{z,t + T)' 



(36) 



(37) 



Since Kt{x,y) in Eq. (32 1 consists only of such time- 
periodic integrals 

Kt{x,y)^Kt+T{x,y). (38) 

Also, since 

dtn{x,t)\t = dtn{x,t)\t+T, (39) 
we conclude from ( [M] ) that 

S{[m,n],x,t) = S{[m,n],x,t + T). (40) 



The first term on the right hand side of Eq. ( 30 ) is also 



periodic with the same period as the density. This implies 
that the entire potential is periodic with the same period 
as the density, i.e.. 



VKsi[m,n],x,t) = VKs{[m,n],x,t + T). 



(41) 



Hence for any possible initial state (labeled by the index 
m) and a density periodic with the period of the external 



field we find that the KS Hamiltonian is also periodic 
with the period of the external field. 

For the Floquet theorem to be applicable in a TDDFT 
framework, the time-dependent KS Hamiltonian must be 
periodic with the period of the external field T = 2tt/uji 
only, i.e., 



VKsi[m,n],x,t) = VKsi[m,n],x,t + T). 



(42) 



Consider now the density being periodic with a period 
T' different from the period of the external field. 



n{x, t) — n{x, t ~\- T'), 



(43) 



as in the above examples in Sees. |IIIB| and [III C[ The 
periods T and T' are incommensurate in general. We 
just have proven that the KS potential is periodic with 
the same period as the density, which implies that 

VKs{[m,n],x,t) = VKs{[m,n],x,t + T'). (44) 

This is in contradiction with the assumption of only one 



period T = 2tt/uji of Eq. (42) which allows the Floquet 
theorem to be applied in the first place. Hence, the Flo- 
quet theorem cannot be applied. 

Here, for our example for which we are able to write 
down an explicit expression for the KS potential, we have 
proven that for any initial KS state it is impossible to 
have a laser-periodic KS potential when the density has 
another period. 



V. CONCLUSIONS 

We investigated the applicability of the Floquet the- 
orem to time-dependent Kohn-Sham Hamiltonians. By 
employing analytically and numerically exactly solvable 
counter examples we showed that, in general, Floquet 
theory is not compatible with time-dependent density 
functional theory. The reason is that, while periodic 
drivers such as laser fields of course render the interact- 
ing many-body Hamiltonian periodic, the corresponding 
Kohn-Sham Hamiltonian, in general, is aperiodic. We 
discussed how the periodicity properties of the single- 
particle density translate to the Kohn-Sham potential. If 
in the Floquet analysis of the many-body time-dependent 
Schrodinger wave function more than one Floquet state 
plays a role — such as for non-adiabatic ramping or res- 
onant interactions — the exact Kohn-Sham potential is 
aperiodic so that the Floquet theorem is inapplicable. 
Further we showed that also the initial-state dependence 
of the time-dependent Kohn-Sham Hamiltonian cannot 
be employed to restore its periodicity. Of course, one 
may view Kohn-Sham-Floquet theory as an approxima- 
tive approach for the study of laser-matter phenomena 
in which resonances and non-adiabaticities are expected 
to be not relevant. 
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